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ABSTRACT 


The  integral  equation 

r 

tP(t}  =  (j.  j  (t-  s}  P(s}  ds,  0  ^  t  <  00  , 

0 

where  jj,  >  0  and  0  <  X  <  1?  has  a  unique  solution  P(t)  which 
satisfies  appropriate  auxiliary  conditions.  Successive 
approximations  to  P(t)  are  derived  by  means  of  a  trapezoidal 
numerical  integration  scheme.  They  converge  uniformly  to  P(t) 
for  0  S  t  <  co  .  Other  approximation  methods  are  described  and 
numerical  examples  are  given.  Applicability  of  the  methods  to 


more  general  integral  equations  is  indicated. 


SOLUTION  OF  AN  INTEGRAL  EQUATION  OF  THE  THIRD  KIND 
BY  SUCCESSIVE  APPROXIMATIONS 

by 

P.  M.  Anselone,  H.  F.  Bueckner  and  D.  Greenspan 

1,  Introduction.  Consider  the  homogeneous  Volterra  integral  equation 

} 

[1.1)  tP(t)  =  (j.  J  (t-s)  P(s)  ds,  0  s  t<  °o  , 

o 

where  p  >  0,  0  <  V  <  1,  and  P[t)  is  a  bounded,  continuous,  non-negative, 
integrable  function  which  satisfies  the  normalization  condition 

00 

[1.2)  /  P(t)  dt  =  I  . 

0 

In  Fredholm's  classification,  [1.1)  is  an  integral  equation  of  the  third  kind. 
Moreover,  the  kernel  has  a  weak  singularity,  Jt  will  be  shown  that  (1.1)  has 
a  unique  solution  P(t)  with  the  prescribed  properties. 

This  paper  is  directed  primarily  to  the  presentation  of  a  method  for  constructing 
successive  approximations  to  P(t)  and  to  proving  that  they  converge  uniformly 
to  P(t)  .  Several  related  approximation  methods  are  indicated  briefly.  The 
methods  are  also  applicable  to  more  general  integral  equations  of  the  form 

t 

[1.3)  tP(t)  =  / K(t~  s}  P(s)  ds,  O^t<oo  . 

0 

However,  treatment  of  such  generalizations  is  deferred  to  another  occasion. 

Sponsored  by  the  Mathematics  Research  Center,  U.  S.  Army,  Madison,  Wisconsin 
under  Contract  No.:  DA-11-022-ORD-2059. 
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Our  interest  in  equation  (1. 1}  was  motivated  by  an  application  [1]  where 
P(t)  is  a  probability  density  function.  This  interpretation  of  P(t)  is  reflected 
in  the  auxiliary  conditions,  which  are  not  independent.  The  following  lemmas 
will  help  to  clarify  the  situation. 

Lemma  1.1.  Let  P(t]  be  a  continuous  solution  of  (1. 1)  .  Then  P(t)  is 
bounded  and  attains  its  maximum  at  some  point  t^  . 

Proof  .  By  (1.1)  , 

t[P(t)|  S  |x 

(1.4)  I P(t)  1 S  max  |p(s)|,  t>  0  . 

OSsgt 

Since  ~  <  1  for  all  t^-^r)1^"  ,  the  lemma  follows, 

l  -  x  1  -  x  ’ 


r  -X 

/(t-s)  ds 

max  [p(s)| 

0 

OSsSt 

Lemma  1. 2.  Let  P(t)  be  a  bounded  solution  of  (1.1)  .  Then  P(t)  is 
integrableon  [0,«>)  and 

(1.5)  P(t)  -  0(t“1_X+€),  t>  0,  0<£S1  . 


Proof .  Suppose  that  [ P (t) {  5  At  *  for  t  >  0,  where  A  >  0  and  0  s  a<  1 
Then,  by  (1. 1)  , 

1 1 P  (t)  |  5  An  /( t-s)  s  ads  , 

0 

|P(t)[  =  A  (j.  B  (1  -  cv  ,  1-X)  t_crX  ,  t>  0  , 
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where  B  is  the  beta  function  .[2]  .  Therefore, 

(1.6}  P(t)  =;  0(t~ff)z^>P(t)  =  0(t~a~V),  0S«<1. 

—OL 

By  hypothesis,  P(t)  =  0(t  )  for  a  =  0  ,  Inductively  let  a=0,  X. ,  2X  ,  ... 

in  (1.6)  to  prove  that  P(t]  =  0(t  m^},  where  m  is  defined  by  (m  -  1)X.  <  1  £  m  X. 

-1+e 

Therefore,  P(t)  =  0(t  }  for  0<  e  £1  .  Another  application  of  (1.6)  yields 

(1.5)  .  So  P(t)  is  integrable  and  the  lemma  is  proved. 

Let  C[0,°o]  denote  the  Banach  space  of  bounded,  continuous,  real  functions 
on  [0,°o]  with  the  supremum  norm.  Let  T  denote  the  integral  operator  defined 
by 

1  r 

(1.7)  (Tf)  (t)  =r/(t-s)  *  f(s)  ds,  0  S  t<  00  , 

0 

for  each  fe  C[0,°o)  such  that  Tf  exists  and  Tf  e  C  [0,oo)  .  Then  (1.1)  can 
be  expressed  as  the  eigenvalue  problem 

(1.8)  TP  =  ~P  . 

|X 

The  operator  T  has  some  interesting  properties  not  commonly  encountered 
in  the  numerical  analysis  of  integral  equations.  Let  f  (t)  be  the  characteristic 
function  of  [e,00)  with  £  >0  in  C1.7)  to  show  that  T  is  unbounded  and, 
hence,  not  completely  continuous.  Since  (1.1)  has  a  non-trivial  solution  for 
every  |x  >  0,  every  positive  number  is  an  eigenvalue  of  T  . 

Since  T  is  a  positive  operator,  it  should  be  expected  that  the  eigenfunction 
P(t)  is  non-negative.  This  is  also  a  direct  consequence  of  the  method  of 
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successive  approximations,  since  every  approximate  solution  is  non-negative. 

Another  interesting  feature  of  the  present  problem  is  that,  as  we  shall  see, 
P(t)  and  all  its  derivatives  vanish  at  t  =  0  .  This  would  seem  to  suggest 
that  only  the  trivial  solution  could  be  obtained  if  the  Volterra  integral  equation 
(1. 1)  were  regarded  somehow  as  an  initial  value  problem.  So  the  usual 
approximations  methods  from  ordinary  differential  equations  are  not  directly 
applicable.  Nevertheless,  the  method  of  successive  approximations  is  based 
to  a  certain  extent  on  such  an  idea  —  each  approximate  solution  vanishes  in 
some  neighborhood  of  t  =  0  . 

The  existence  of  P(t}  will  follow  from  the  method  of  successive 
approximations.  The  uniqueness  is  established  in  Section  2  by  means  of  Laplace 
transforms.  The  Laplace  transform  approach  can  be  used  also  for  the  existence 
of  a  solution,  but  it  is  difficult  to  show  by  this  means  that  the  solution  is  non¬ 
negative  or  has  other  significant  properties.  Moreover,  the  method  of  successive 
approximations  is  applicable  to  more  general  equations  (1.  3)  for  which  the  Laplace 
transform  of  P(t)  may  not  be  explicitly  derivable. 
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2.  Laplace  Transform  Analysis..  Assume  that  P(t)  is  a  solution  of  (1.1) 
with  the  prescribed  properties.  Then  the  Laplace  transform  of  P(t)  , 

00 

(2.1)  P(z)  =  /  e  Zt  P(t)  dt  (z  =  x  +  iy)  , 

0 

is  defined  for  all  xgO  .  Transform  (1.1)  and  use  the  convolution  theorem  to 
obtain 


(2.2)  P'(z)  =  -(ir(l-X)zV“'1  P(z)  ,  x>  0  , 

with  the  principal  branch  of  zX  1  .  By  (1. 2)  , 

(2.3)  P(0)  =  1  . 


The  unique  solution  of  (2.2)  and  (2.3)  is 

X 


(2.4) 


,  ~«z 

P(z)  =  e 


n 


urg-x) 


with  the  principal,  branch  of  z  ,  Since  P(z)  determines  P(t)  uniquely  (e.  g.  , 
by  means  of  the  complex  inversion  integral),  P(t)  is  determined  uniquely  by 


(1.1)  and  (1.2)  . 

When  convenient,  P(t;  X  ,  p.)  will  be  written  for  P(t)  to  indicate  the 
dependence  on  X  and  p.  .  We  assert  that 


(2.5)  P(t;  X;fx)  =(^}1/X  P(£)1/X  t;X,  v)  . 

H*  r 

One  way  to  prove  this  is  to  show  that  both  members  of  (2.5)  satisfy  (1.1)  and 
(1.  2)  .  Another  way  is  to  show  that  both  members  have  the  same  Laplace  transform, 
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given  by  (2.  4)  .  In  view  of  (2.5)  ,  it  suffices  to  determine  P(t;  X  ,  (i  )  for 
any  single  fixed  value  of  (j.  . 

In  particular,  (2.4)  yields 


(2.6) 


p(z)  =  e-^  if  F  =  i/r|i-i]  . 


Standard  works  on  Laplace  transforms  [3],  [4]  give  the  inverse  transforms  of 

-z^  1  2 

e  for  X.  =  —  and  X  =~  .  Thus, 


(2.7) 


P(t; 


1 


2  ; 


2*/  ir 


:) 


-l/4t 


and 


(2.8) 


2 

3  r  (1/3) 


)  = 


-2/27t2 

e 

2^  3 it  t 


W 


-1/2,  -1/6 


(-4/27t2) 


where  W  denotes  a  Whittaker  function.  These  results  and  (2.5)  yield 
P(t;  1/2,  (j.)  and  P (t;  2/3,  jx)  for  all  (j.>  0  . 

Pollard  [5]  derived  a  series  expansion  for  the  inverse  transform  of 
which  yields 


> 


(2.9)  P(t;  V,  X/r(l-X))  =  Z 

n=l 


n+1 


sinfnTrX  )  r  (nX  +  1) 


trn 


,  nX  + 1 


Again,  (2.5)  yields  P(t;  X,  (j.)  for  all  X  and  p.  .  This  series  is  useful  for 
calculating  P(t)  for  large  values  of  t  .  However,  it  converges  very  slowly  if 
t  is  only  of  moderate  size. 
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Van  der  Corput  [6]  studied  the  asymptotic  behavior  as  t-*-0  of  the  sum  of 
the  series  in  (2.9)  .  His  formulas  are  useful  for  evaluating  P(t)  for  small  t  , 
Moreover,  they  imply  that  P(t)  and  all  its  derivatives  vanish  at  t  =  0,  as  was 
mentioned  earlier.  Note  that  (2.7)  and  (2.8)  are  consistent  with  this  result. 

Another  method  of  calculating  P(t)  is  by  means  of  the  complex  integral 
formula  for  the  inverse  Laplace  transform.  Thus,  by  (2.4), 

§  +  iig  % 

(2.10)  P(t)=~lim  /  eztQ-nz  dz,  £>0  , 

rj-»oo  |-ir| 

where 


(2.11) 


X  V  iX.0 
z  =  r  e 


for  z  =  re' 


i0 


TT  TT 

0<r 


By  the  reflection  principle, 


(2.12)  P(t)=^/  +  1  Re(eZt"  ** 

e+io  [ 


which  is  an  improper  Riemann  integral.  Let  z  =  %  +  iy  to  obtain 


|  dz,  £  >  0  , 


00  \ 

(2.13)  P(t  )-~f  nr  cos  cos  [yt  -  £2r^  sin  \  0]  dy,  §>0  , 


where 


, ,2  2,1/2 

r  =  (l  +  y  }  , 


°  =  0<I 


(2.14) 


0  =  arctan  [y/£], 
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Let  §-*-0  to  obtain 

00  X 

(2.15)  P(t)  -  J  e  ^  cos(^Tl/2')  cog  _  Qy*-  sinCXir/Z)]  dy 
17  0 

Although  P(t)  can  be  approximated  by  integrating  (2.13)  or  (2.15) 
numerically,  a  rather  long  calculation  is  involved  for  each  value  of  t  .  The 
method  of  successive  approximations  presented  below  provides  an  alternative  to 
the  use  of  these  integrals  or  the  formulas  of  Pollard  and  Van  der  Corput. 
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3.  The  Function  R(t)  .  The  method  of  successive  approximations  will 
be  presented  first  for  a  solution  of  equation  (1.1)  which  satisfies  a  different 
normalization  condition.  Thus,  consider  the  integral  equation 

(3.1)  t R(t)  =  (j.  J  (t-s)  R(s)  ds,  0£  t<  w  , 

0 

where  R(t')  is  bounded,  continuous,  non-negative,  integrable,  and 

(3.2)  max  R(t)  =  1  . 

0gt<°3 

(The  maximum  is  attained  by  Lemma  1.1  .)  The  functions  P(t)  and  R(t)  are  related 


by 

00 

(3.3) 

P(t)  =  R(t) /  f  R'(s)  ds  , 

0 

0  S  t  <  °°  , 

(3.4) 

R(t)  =  P(t)/max  P(s)  , 

0  £s<°° 

0  S  t  <  «3  r 

The  existence  of  both  P(t)  and  R(t)  will  follow  from  the  method  of 
successive  approximations.  The  Laplace  transform  analysis  yields  the  uniqueness. 


I 
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4.  The  Method  of  Successive  Approximations.  For  h  >  0,  let  Rk(t)  >■ 

0  £  t  <  °° ,  denote  a  continuous,  non-negative,  bounded,  piecewise-linear 
function  with  possible  changes  in  slope  only  at  the  points  t  =  nh,  n  =  1,,  2,  . , 
Thus,.  R^Ct),  t  50,  is  determined  in  terms  of  the  values  R^(nh},  n^0,  by 

(4.1)  Rh(t)  =  tnty  ~t'Rh(nh)  +  Rh(Cn+l)h),  nhgtS(n+l)h,  n  go  . 


Assume  that  R^t)  satisfies  equation  (3.1)  at  the  points  t  =  nh: 


(4.2) 


nh 

nh  (nh)  -  [i  J  (nh-s)  R^fs)  n  =  °  • 


Other  conditions  will  be  imposed  later.  Ultimately,  we  shall  let  h-*-0  . 

The  procedure  for  deriving  R^(t)  may  be  called  a  trapezoidal  method  by 
analogy  with  the  trapezoidal  rule  for  numerical  integration.  Other  methods,  which 
correspond  to  other  numerical  integration  rules,  will  be  discussed  in  Section  7. 

By  (4. 1)  and  (4. 2)  , 

n-1  (k+l)h 

nh  R,  (nh)  =  p.  £  /  (nh-s)  R,(s)  ds  , 
k=0  kh 

-X.  f  n-1  r 


nRhCnh)  =(1-M(2-X)  \  k^Q 


(2  -  X  )(n  -  k)1_X  -  (n-k)2"V  +  (n_k_1}2-\ 


Rh(kh) 


n-1 

+  E 

k=0 


(n  -  k)2"X  ~{n  -  k-l)2_X  -(2-X  )(n-k-l)1_X 


Rh((k+l)h) 


Replace  k  by  k  - 1  in  the  second  sum  and  collect  terms  to  obtain 


(4.3)  (n-nh)  Rh(nh)  =  nh 


rW°i+%  VkVkh> 

k=  0 


ngl, 


where 

\ 

(4.4)  bn  =  (2-X)  n1_X  +  n2"X  -  (n  +  l)2^  ,  ngl, 

(4.5)  cn  =  (n  +  1)2"K  -  2n2"X  +  (n  -  1)2"X  ,  nfcl, 


and 


(4.6) 


u.h  ^ 

nh  =  (1-X)(2-X) 


For  future  reference,  note  that 

(4.7)  (l-nh)  Rh(h)  =  nh(l-\)  Rh(0)  , 

(4.8)  nh  T 00  as  h  4r  0  , 
and 

(4.9)  cn  >  0  ,  nil, 

2-X 

since  each  cn  is  a  second  difference  of  a  function  t  with  a  positive  second 
derivative. 

If  n^  is  not  an  integer,  then  (4.  3)  determines  R^(nh),  ngl,  inductively 
in  terms  of  R^(0)  .  Then  (4.1)  yields  a  solution  (t)  of  (4.2)  .  If  n^>l 
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hypothesis  ,  it  follows  that  R^(0)  and,  hence,  R^(t)  =  0  if  h  is  sufficiently 
small.  This  case  will  not  be  considered  further. 

We  assume,  henceforth,  that  n^  is  a  positive  integer,  i.  e. ,  we  restrict 
h  to  the  bounded  and  countable  set 

(4.10)  H  =  {h:  n  =  1,  2,  3,  ...  >  . 

Nrw  (4.3)  determines  R^(nh)  for  IS  n<  n^,  but  not  for  n^,  in  terms 
of  •  As  before,  (4.7),  (4.8)  and  R^(t)S0  imply  that  R^(0)  =  0  if 

n^>l  .  Let  Rja(°)  =0  for  all  he  H  .  Then ,  by  (4.3)  , 

(4.11)  R^(nh)  =  0,  0Sn<  n^ 

Assume  that 

(4.12)  Rh(nhh}>0  . 

Since  (4.  3)  is  now  satisfied  automatically  for  n  s  n^,  it  reduces  to 

n-1 

(4. 13.)  (n  -  nh)  R^nh)  =  nh  £  Cn_k  R^kh),  n  >  nh  , 

k=n, 

h 

which  determines  R^(nh),  n  >  n^,  inductively  in  terms  of  R^n^h)  .  Then 
(4.1)  yields  a  solution  R^t)  °f  (4.2)  .  Since  (4.1)  and  (4.13)  are  linear 
relations , 

(4.14)  Rh(t)  =  Rh(nhh)  ^(t)  > 

where  R^(t)  is  the  particular  solution  with  R^(n^h)  =  1  . 
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Equations  (4.1)  and  (4.11)  yield 

(4.15)  Rh(t)  =  0,  0|t£(nh-  1)  h  , 

where  (n^-l)h-*0  as  h-*-0  by  (4.6)  .  Thus,  both  R^(t)  and  R(t)  are 
"very  flat"  at  t  =  0  .By  (4. 1),  (4.  9),  (4. 12)  and  (4. 13), 

(4.16)  Rh(t}>0,  t>(nh-l)h  . 

Therefore,  the  hypothesis  that  (t)  SO  is  satisfied. 

Lemma  4. 1.  R  (t)  is  bounded  and  attains  its  maximum  at  some  point 
t  =  khS(i:^-)lA  . 

Proof.  The  proof  is  based  on  (4.1)  and  (4.  2)  .  It  is  analogous  to  that 
for  Lemma  1.1. 

Thus  far,  R^(t)  is  the  general  non-trivial  solution  of  (4.2)  which  satisfies 
the  given  conditions.  By  analogy  with  (3.  2)  assume  now  that 

(4.17)  max  R  (t)  =  1  . 

0^t<  oo 

Then  R^(t)  is  determined  completely  and  is  given  explicitly  by 

R^Ct)  =  R^Ct)/  max  r!(s)  , 

0  Ss<  oo  . 


(4.18) 
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This  formula  is  convenient  for  calculation.  In  view  of  Lemma  4. 1,  the  maximum 
of  R^(t)  is  a  computable  quantity.  Each  function  (t) ,  h  e  H,  is  a 
continuous,  piecewise-linear,  non-negative  approximate  solution  of  (3.1}  and 
(3.2)  . 

Lemma  4.  2.  R^(t)  is  integrable  on  [0,eo)  .  Moreover,  for  0  <  e  £  1  , 

there  exists  A(«)  >  0,  independent  of  h,  such  that 

(4.19)  IL(t)£  A(€)t"1_X+€,  t>0,  h  c  H  . 

Proof.  The  proof  is  based  on  (4.1)  and  (4.2)  .  It  is  analogous  to  that 
for  Lemma  1.  2. 

Since  R^(t)  is  piecewise  linear  and  11^(0}  -  0, 

00  00 

(4.  20}  /  R^t)  dt  =  h  Yj  Rj-Cnh)  • 

0  n=0 

In  view  of  (4.19),  this  quantity  can  be  calculated  to  any  desired  accuracy. 
Define 

CO 

(4.21)  Ph(t)  =  Rh(t}/ /  Rh(s)  ds,  t  SO  . 

Then 

nh 

(4.22)  nh  P^nh)  =  (j ,  J  (nh-s)  P^(s)  ds,  nSO  , 


and 


Each  function  P^(t),  h  «  H,  is  a  continuous,  piecewise-linear  non-negative 
approximate  solution  of  (1. 1}  and  (1. 2}  . 
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5.  Properties  of  the  Approximations. 

Lemma  5.1.  The  functions  (t) ,  h  e  H,  are  uniformly  equicontinuous  on 
each  interval  t^  £  t  <  °°  with  t^  >  0  .  Moreover, 

(5.1)  Iiyt)  -Rh(s)|  g  J 


?2\+l 

— (t-  s)  +  (t-s}>  ,  t  Ss>  0,  he  H 


Proof .  For  ngmgO,  it  follows  from  (4.2)  and  (4.17)  that 

nh  mh 

[nhR^tnh)  -  mhR^tmh)!  =  |j.  1  J  (nh-s)  R^(s)ds  -  J  (mh-s)  ^  R^(s)  ds| 

{  nh  -v  mh  -\ 

=  (j.  <  /  (nh-s)  ds  +  J  ((mh-s)  -(nh-s)  ]  ds 

mh  0 


-  _i±_ 
"  1-X. 


2(nh-mh)1“XL  +  (mh)1^  -  (nh)1"V 


^  ^-(nh-mh)1_V 


By  symmetry, 


(5.2)  [nhR^(nh)  -  mhR^mh}!  §  [nh  -  mh|3'  K  ,  m,  nSO 
It  follows  from  (4.17)  and  (5.2)  that 

[  nh  R^(  nh)  -  nh  Rh(  mh)  [  s  [  nh  R^(  nh}  -  mh  R^(  mh)  [  +  [  mh  R^  (  mh)  -  nh  R^(  mh )  [  , 

(5.3)  lRh(nh)  -  ^(mh)!  £  ^  |  [nh-mh[1X'  +  |nh-mh|j>  ,  n>0  . 


t 


I 
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Let  nhg  s,  t£(n+l]h  .  Then,  by  (4.1)  and  (5.3), 


i  v*>  -  Rh'sn  =  ^ 


V 

It  -  s[ 


[Rh((n  +  l)h)  -  R^(nh)| 


1 


u1_X- 


h  (n4  l)h  1  -  X 


h*  "  4  h  l  = 


JjL 


(n+l)h  1  1  -  V  hX 


4  It-; 


SO  that 


(5.4) 


I R^(t)  -  Rh(s)[  £  ^  |t-  si1  X  4  It-  s[  |  ,  nh  g  s,t  g  (n4l)h,  s>  0 


Therefore,  (5.1)  is  valid  in  this  case. 

Now  assume  that  s  g  nh  g  t  for  some  n  .  Define  m  and  n  such  that 
(m  -  l)h  g  s  g  mh  g  nh  g  t  g  (n+  l)h  .  By  (5.3)  and  (5.4)  , 


\{t)  -  Rh(s)[  g  \(t)  -  Rh(nh)|  +  [^(nh)  -  Rh(mh)|  4  ^(mh)  -  Rh(s}| 


=  ~  |  j~^-(t  -  nh)1  V  4  (t  -  nh)  i  4  ”  j  ^^-(nh  “  mh^  V  +  (nh  -  mh) 


+  s  I  A"  ^mh  "  s)2  ^  +  (mh  “  s)  |  > 


(5.5)  I\(t)-Rh(S)i  =  s  Irt 


,1-X 


1-X 


(t-nh)  4(nh-mh)  4(mh-s) 


,1-X 


4  (t  -  s) 


Since 


0<\<  1,  the  function  f(t)  =  t1  ^ ,  t§0,  is 


1-X  ,  1-X  /  , 

a  4b  <  ajhjb 


1-X 


1-X  ,  1-X 

a  4  b 


concave: 

g  ZX  (a4b)1_V 
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for  a,  b  SO  .  It  follows  that 


(5.  6)  a1-X  +  b1_X  +  c1_X  ^  22X  (a  +  b  +  c)1  X  ,  a,  b,  c  2.0  . 

Then  (5.5)  and  (5.6)  yield  (5.1)  for  this  case.  Therefore,  (5.1)  is  valid  in 
general  and  the  lemma  is  proved. 

Several  properties  of  the  constants  c^  defined  by  (4.  5)  will  be  needed  in 
the  proof  of  the  next  lemma.  First,  note  that 


c  =  n 
n 


2-\ 


[<1  +  n 


>2"X  -  2  +  (l 


1  2-V 

Use  three  terms  of  the  binomial  expansions  for  (1+  — )  to  prove  that 


(5.7) 


nXcn->(l-\)(2-V) 


as  n  • 


oo 


Next,  note  that  c^  =  c(n),  where 

c(t)  =  (t+l)2“X  -  2tZ_X  +  (t-l)2"X,  tSl  . 


Then 

~^-c‘(t)  =  (t+l)1_X  -  2t1_X  +  (t  -  1)1“X  ,  til  . 

Since  this  is  a  second  difference  of  the  function  g(t)  =  t^  X  ,  t  SI,  which  has 

a  negative  second  derivative,  c'(t)  <  0  .  It  follows  that  both  c(t)  and  {c  } 

n 

are  monotone  decreasing.  Therefore,  by  (5.7)  , 

c  J,  0  as 
n 


(5.8) 


n  -*•  co  . 
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Lemma  5,2.  Let  0  <  .  Then,  for  all  sufficiently  small  he  H  , 


the  functions  R,  (t)  are  monotone  non-decreasing  on  the  interval  0£  t  2  t 


1  • 


Proof.  In  view  of  (4.11),  (4.12)  and  Lemma  4.1,  there  exists  a  unique 
integer  for  each  h  «  H  such  that 


(5.9) 


Vkh)  <  \mi)h),  nh  -  1<  k  <  Nh  , 


Rh(Nhh)  Vh  +  1W 


Then  R^(t)  is  non-decreasing  for  0  5  t  S  h  .  By  (4.13)  , 


( Nh  + 1  -  nh)  Rh  C(Nh  +  l)h)  -  (Nh  -  nh)  Rh(Nh  h) 


V  iVh  v^'h 

N, 


N, -1 

h 


-  nH  E  CTKT  4-1-lf  Rh^kh^  "  nh  E  CN  _k  Rh^kh^  • 

k=0  h 


Ti  '"'N.+l-k  V 
k=0  h 


Replace  k  by  k  +  1  in  the  first  sum  and  use  R,(0)  =  0  to  obtain 


(Nh+l-nh)[Rh((Nh+l)h)  -  Rh(Nhh)]  +  Rh(Nhh) 


Nh"1 

=  nh  l  c  k  [Rh((k+l)h}  -  Rh(kh)]  . 


Hence,  by  (5.8),  (5.9)  and  Rk(0}  =  0  , 


Nh-‘ 


W,Snh°N  Z  |Il(Wl.-V««-1,XW 
h  k=0  h 
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Therefore,  .1^  c  £  1  .  By  (4.6)  , 
h 


(5.10) 


(1-MU-M  h\ 


so  that  -*-0  as  h  0  .  In  view  of  (.5.8),  -*•  00  as  h-*-0  .  By 

h 

(5.10)  and  (5.7)  , 


^  °N. 

(Nh  h)  a  H  » 


h-0  , 

N.  -►  00  . 


Therefore,  h  2t^  for  h  sufficiently  small.  Since  Rh(t)  is  non-decreasing 
for  0,S  t  £  h,  the  lemma  is  proved. 

We  conclude  this  section  with  the  statement  of  an  auxiliary  result  of  a 
general  nature  which  will  be  used  shortly. 

Auxiliary  lemma  .  If  a  sequence  of  monotone  functions  converges  pointwise 
to  a  continuous  function  on  a  closed  interval,  then  the  convergence  is  uniform. 

This  is  an  old  result  whose  proof  is  not  difficult.  For  a  recent  reference, 


see  [7], 
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6.  The  Convergence  of  the  Successive  Approximations  .  We  consider  first 
the  convergence  of  the  functions  (t)  as  h 0  through  H  .  We  shall  regard 
{R^:  h  *  H}  as  a  sequence  ordered  by  letting  h  decrease  through  H  or, 
equivalently,  by  letting  n^  increase  through  the  positive  integers. 


Lemma  6. 1.  Every  subsequence  of  {R^:  he  H)  has  a  further  subsequence 
which  converges  uniformly  on  each  finite  interval. 


Proof.  Since  0  <_  R^(t)_<_  1,  Lemma  5. 1  implies  that  the  Arzela-Ascoli 
theorem  is  applicable,  and  by  that  theorem  there  exist  successive  (nested) 
subsequences  of  (R^:  h  e  H}  which, respectively,  converge  uniformly  on  the 
intervals  ,  m];  m  =  E,  3,  4,  ...  .  Then  the  usual  diagonal  procedure  yields 
a  single  subsequence  which  converges  uniformly  on  every  one  of  the  intervals 
[~,  m]  .  Since  R^(0)  =  0  for  all  he  H,  the  subsequence  converges  pointwise 
for  t  >  0  .  It  follows  from  Lemma  5.  E  and  the  Auxiliary  Lemma  that  the  subsequence 
converges  uniformly  on  each  finite  interval. 


Lemma  6.  E  .  Every  subsequence  of  (R^:  h  e  H}  which  converges  uniformly 
on  each  finite  interval  converges  to  a  solution  R(t)  of  (3.1)  with  the  prescribed 
properties. 


Proof.  Suppose  that  ( t)  -*■  R(t)  as  h  -*  0  through  a  subset  H*  of  H  . 

The  assumed  uniform  convergence  and  the  continuity  of  each  function  R^(t)  imply 
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that  R(t)  is  continuous.  Since  0  £R^(t)  si,  we  have  0  g  R(t)  g  1  .  By 
^4.17)  and  Lemma  4.1  , 

max  R(t)  =■  lim  max  R,  (t)  =  1  , 

l-o  ^ 

-^•1-X  j  h«  H*  1  X 

so  that  max  R(t)  =  1  .  In  (4.2),  let  n->°° ,  h  -*■  0  and  nh  -*•  t  with  h  «  H1 
0gt<oo 

to  obtain  (3.1).  Finally,  R(t)  is  integrable  by  Lemma  1.  2  .  Thus,  the  lemma 
is  proved. 

Theorem  6.1.  There  exists  a  unique  solution  R(t)  of  (3.1)  with  the 
prescribed  properties,  and  R^(t)  -»  R(t)  uniformly  for  OSt<oo  as  h->0  through 

H  . 


Proof.  The  existence  of  R(t)  follows  from  Lemmas  6.1  and  6.2.  The 
uniqueness  is  a  consequence  of  (3.4)  and  the  uniqueness  of  P(t),  which  was 
established  in  Section  2.  Fix  t  e  [0,°°)  arbitrarily  and  let  r(t)  be  any  limit 
point  of  the  numerical  sequence  (R^(t):  h  e  H}  .  Then  the  sequence  of 
functions  (R^:  h  «  H)  has  a  subsequence  which  converges  at  t  to  the  value 
r(t)  .  By  Lemmas  6.1  and  6.2,  there  is  a  further  subsequence  which  converges 
to  R  on  [0, oo)  .  Therefore,  r(t)  =  R(t)  and,  hence,  R^(t)  -*■  R(t)  as  h-*-0 
through  PI  . 


I 
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The  pointwise  convergence  of  R^(t)  to  R(t)  and  the  uniform  equicontinuity 
established  in  Lemma  5.1  imply  that  R^(t)  -*■  R(t)  uniformly  on  each  interval 
[t  ,  t  ]C(0,°o)  .  By  Lemma  5.2  and  the  Auxiliary  Lemma,  the  convergence  is 
uniform  on  each  finite  interval.  Finally,  the  convergence  is  uniform  for 
0  S  t  <  «  by  (4. 19)  . 

Theorem  6,2.  There  exists  a  unique  solution  P(t)  of  (1.1)  with  the 
prescribed  properties,  and  P^(t)*^P(t)  uniformly  for  0  §  t  as  h-*-0 
through  H  . 

Proof.  The  existence  of  P(t)  follows  from  (3.3)  and  the  existence  of 
R(t)  .  The  uniqueness  was  established  in  Section  2.  By  Lemma  4.2,  Theorem 
6.1,  and  the  Lebesgue  dominated  convergence  theorem 
oo  oo 

(6.1)  f  R,  (t)  dt -*■  J  R(t)dt  as  h  0  . 

0  n  0 

Therefore,  by  (3.  3),  (4.  21)  and  Theorem  6.1,  P^(t)  converges  uniformly  to 
P(t)  and  the  theorem  is  proved. 

In  addition,  it  can  be  shown  that  P^(t)->P(t)  in  L^(0,  oo)  as  h->-0 
through  H  .  This  is  important  if  P(t)  is  interpreted  as  a  probability  density. 
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7.  Related  Methods  for  Deriving  Successive  Approximations  .  Though  our 
remarks  in  this  section  are,  for  the  sake  of  brevity,  restricted  primarily  to 
R^(t),  analogous  results  hold  for  P^(t)  • 

The  trapezoidal  method  described  in  Section  4  has  the  following  features. 

(a)  R^(t)  is  determined  by  linear  interpolation  from  the  values 
H^Cnh),  n  >  0  ; 

(b)  \{t)  satisfies  the  integral  equation  (3.1}  at  the  points 
t  =  nh,  n  >  0  ; 

(c}  R,  (t)  satisfies  the  normalization  condition  (3.2)  . 

Variants  of  the  method  are  obtained  if  other  interpolation  rules  are  used 
in  (a)  .  For  example,  consider  R^(t)  a  step  function  such  that 


(7.1)  R^(t)  =  R^(nh),  nh£t<(n  +  l)h,  ngO  . 

It  would  follow  from  (b)  that 


(7.2)  nRh(nh) 


K=U 


n  gl 


Then  (7.2)  determines  R^(nh),  ngl,  by  induction  in  terms  of  R^(0),  R^(t) 
is  obtained  from  (7.1),  and  R^(0)  is  chosen  so  that  (4.17)  is  satisfied. 
Theorem  6. 1  is  valid  also  for  this  case.  Note  that  R^(0}^0,  in  contrast  with 
the  trapezoidal  method. 

As  another  example,  consider  R,  (t)  a  step  function  such  that 


(7.3) 


=  ^((n  +  i-Jh),  nh<t  §  (n  +  1)  h, 


ngO  . 
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Then  (b}  yields,  for  each  nil  , 


(7.4) 


X.  n-1 


(n-‘TT)Rhl"hi  =  ftr  2  Vkh> 

k=l 


(n-k  +  l)1^  -  (n-k)1-X 


Here,  analogous  to  the  trapezoidal  method,  we  assume  that  m^  =  is  an 

integer,  and  let  R^(nh)  =  0  for  Oin<  m^  .  The  rest  of  the  derivation  of 
R^(t)  is  similar.  Again,  Theorem  6. 1  is  valid. 

A  third  and  last  example  involves  Rh^  determined  by  parabolic  interpolation 
from  three  successive  values  R^(nh)  .  Once  again,  Theorem  6. 1  is  valid. 

Further  details  are  omitted. 

Finally,  it  is  worth  noting  that  the  original  equation  (1.1}  provides  a 
natural  means  of  smoothing  an  approximate  solution  R^Ct)  • 


(7.5) 


i> 


The  convergence  of  R^£t)  as  h-*-  0  is  not  investigated  in  this  report. 


-26- 


#345 


8.  Examples.  The  numerical  calculation  performed  in  support  of  the  methods 
described  in  this  paper  were  too  extensive  for  complete  reproduction  and  though 
a  variety  of  values  of  X  and  p.  were  considered,  only  two  typical  examples  in 
which  X  =  ~  ,  p.  =  — - —  will  be  discussed.  Programming  was  performed  by 

2  2  V7 

J.  A1  Abdulla  and  calculations  were  completed  on  the  CDC  1604. 


Example  1.  For  X  =  p.  ,  the  exact  solution  of  {1.1}  -  {1.2}  is  given 

2 's/  it 

by  (2.7)  and  (2.13)  .  The  expression  in  (2.13}  was  approximated  by 


(8.1) 


1  f  j 

—  I  exp  <  t - 
W  0  1 


^  2,1/4 

(i  +  y  ) 


,  ,,  2, -1/2 
i+ii+ij.,,,. 


cos  i  yt-(l  +  y2}^4 


irXi.tx2)."1/2 


1/2 


dy 


Then,  (8.1)  was  evaluated  for  a  variety  of  values  of  t  in  the  range  OS  t  §  1 
by  means  of  Simpson's  rule  with  h  =  .  0001  .  It  was  found  that  any  (3  S4500 
yielded  approximations  which  agreed  with  the  exact  solution  to  at  least  five 
decimal  places. 

Example  2.  For  X  =  7,  p.  =  “ — I —  ,  the  exact  solution  of  (3.2)  -  (3.  3)  is 

L  2  VtT 

(8.2)  R(t)  =  P (t) /  max  P(s)  , 

0Ss<°o 

where  P(t)  is  given  by  (2.7)  .  For  n^  =  6,  h  is  approximately  .0039  by 
(4.6),  and  R^(t)  and  R^(t)  were  generated  by  means  of  (4.13)  and  (4.17)  . 
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Selected  results  are  displayed  in  Table  I.  It  must  be  emphasized  here  that  an 
unexpected  stability  always  accompanied  the  trapezoidal  method.  In  all  our 
examples  roundoff  error  never  accumulated  excessively  and  R^(t}  always 
appeared  to  damp  out  even  though  n  was  taken  to  be  as  large  as  100,  000. 
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Table  I 


n 

t  =  nh 

Rh<‘> 

v*> 

R(t) 

6 

.  0234 

1.0000 

.  0012 

.00195 

7 

.  0273 

4.9706 

.  0064 

.  00713 

10 

.  0390 

51.6291 

.  0670 

.  06512 

14 

.  0546 

193.8953 

.  2516 

. 24547 

20 

.  0780 

443.4639 

.  5755 

. 56781 

25 

.  0975 

599.  1857 

.  7776 

.77131 

30 

.  1170 

694.4550 

.9012 

.89961 

38 

.  1482 

763. 3806 

.9907 

.98952 

42 

.  1638 

770.  5  151 

1.  0000 

1.  00000 

45 

.  1755 

768. 5804 

.9974 

.99732 

50 

.  1950 

755.  8736 

.9810 

.98282 

54 

.2106 

740. 0231 

.9604 

.96290 

60 

.  2340 

710.8442 

.9225 

.92576 

62 

.  2418 

700. 2745 

.9088 

.91223 

65 

.  2535 

683.9963 

.  8877 

.  89  135 

68 

.  2652 

667.4421 

.  8662 

.  87006 

70 

.  2730 

656. 3554 

.8518 

. 85578 

78 

.  3042 

612. 5302 

.7949  . 

.79921 

86 

.  3354 

570.7857 

.  7407 

.  745  16 

94 

.  3666 

531.9508 

.  6903 

. 69480 

102 

.  3978 

496. 2682 

.  6440 

. 64846 
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